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$_i ■ the collapse of subcritical bubbles in models with symmetric and asymmet- 

ric double-well potentials. By a combination of analytical and numerical 
work we explain several of their properties, including the conditions for 
their existence, their longevity, and their final demise. We discuss several 
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nucleation during cosmological phase transitions may have wide-ranging 
consequences. 
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I. INTRODUCTION 



The search for static, localized, non-singular solutions of nonlinear field theories has 
by now a long history [[L]]. In (l+l)-dimensions, it is possible to find exact static solutions 
to the nonlinear Klein-Gordon field equations for certain interacting potentials, such as 
the kink solutions of sine-Gordon or 4 models. For a larger number of spatial dimensions, 
Derrick's theorem forbids the existence of static solutions for models involving only real 
scalar fields There are several ways to circumvent Derrick's theorem, by invoking 
more complicated models with two or more interacting fields. Well-known examples 
include topological defects such as the 't Hooft-Polyakov monopole or the Nielsen- Olesen 
vortices 0. Topological conservation laws guarantee the stability of these configurations. 

It is also possible to find localized time-dependent but non-dissipative solutions of 
nontopological nature, the so-called nontopological solitons The simplest model of 
a nontopological soliton in the context of renormalizable theories has a complex scalar 
field quadratically coupled to a real scalar field with quartic potential. The stability of 
the configuration comes from the conserved global charge Q carried by the complex field 
which is confined within a spherically-symmetric domain formed by the real scalar field. 
One can show that for Q larger than a critical value, the energy of the configuration 
is smaller than the energy of Q free particles. There has been a recent upsurge of 
interest on nontopological solitons due to their potential relevance to cosmology and 
astrophysics || . If one waives the requirement of renormalizability, it is possible to find 
nontopological solitons for models with a single complex scalar field, by invoking, e.g., a 
(f) 6 term in the potential. These are the so-called Q-ball solutions discovered by Coleman 
and collaborators ||. 

In the present work we will go back to the simple models involving only a self- 
interacting real scalar field and study the properties of time- dependent spherically- 
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symmetric solutions. Due to the constraint imposed by Derrick's theorem, these con- 
figurations have been somewhat overlooked in the literature (but not completely, as we 
will discuss below). Why should anyone bother with solutions which are known to be 
unstable? One possible answer is that instability is a relative concept, which only makes 
sense in context, that is, when the lifetime of a given configuration is compared with 
typical time-scales of the system under study. Thus, unstable but long-lived configu- 
rations may be relevant for systems with short dynamical time-scales. Another answer 
is that a detailed study of these configurations can greatly clarify dynamical aspects of 
nonlinearities in field theories and the role they play in several phenomena, ranging from 
nonlinear optics to phase transitions both in the laboratory and in cosmology J7|. 

One of the motivations for studying the evolution of unstable spherically-symmetric 
configurations comes from the work of Gleiser, Kolb, and Watkins on the role subcritical 
bubbles may play in the dynamics of weak first order phase transitions pi. Considering 
models with double-well potentials in which the system starts localized in one minimum, 
these authors proposed that for sufficiently weak transitions correlation-volume bubbles 
of the other phase could be thermally nucleated, promoting an effective phase mixing 
between the two available phases even before the critical temperature is reached from 
above. This could have important consequences for models of electroweak baryogenesis 
which rely on the usual homogeneous nucleation mechanism |J. However, Gleiser, Kolb, 
and Watkins did not include the shrinking of the bubbles in their estimate of the fraction 
of the volume occupied by each of the two phases, leading some authors to question 
their results [[Hj. Since then, Gleiser and Gelmini included the shrinking of the bubbles 
into the original estimates, concluding that for sufficiently weak transitions subcritical 
bubbles are indeed nucleated at a fast enough rate to cause substantial phase mixing 



Although an improvement, the modeling used to describe the bubble shrinking was 



still too simplistic, as it assumed that the bubbles just shrunk with constant velocity. 



The evolution of spherically-symmetric unstable solutions of the nonlinear Klein- 
Gordon equation was originally studied numerically in the mid-seventies by Bogolubsky 
and Makhankov Using a quasiplanar initial configuration for the bubbles (that is, 
a tanh(r — Rq) profile, with Rq the initial radius), these authors discovered that for a 
certain range of initial radii the bubble evolution could be described in three stages; 
after radiating most of its initial energy the bubble settled into a regime which was 
quite long-lived, with a lifetime which depended on the initial radius. The bubble then 
disappeared by quickly radiating away its remaining energy. These configurations were 
called "Pulsons" by these authors, due to the pulsating mechanism by which they claimed 
the initial energy was being radiated away. Their results were recently rediscovered and 
refined by one of us After a more detailed analysis of these configurations, it 

became clear that their most striking feature was not the pulsating mechanism by which 
bubbles radiate their initial energy, but the rapid oscillations of the field's amplitude at 
the core of the configuration during the pseudo-stable regime, in a manner somewhat 



analogous to resonant breathers in kink-antikink scattering [II]. In fact, it was realized 
that during the pseudo-stable regime almost no energy is radiated away and the radial 
pulsation is actually quite small in amplitude. Hence the name "Oscillon" was proposed 
instead. It was also shown that these configurations appear both in symmetric and 
asymmetric potentials, are stable against small radial perturbations, and have lifetimes 
far exceeding naive expectations. However, not much else has been done in order to 
explore the properties of these configurations. Other works on this topic were concerned 
in establishing the existence of these solutions for other potentials, such as the sine- 
Gordon and logarithmic potentials, different symmetries, and somewhat limited stability 
studies [|T5|j . 

By a combination of analytical and numerical methods, we will shed some light on the 
properties of these configurations (henceforth oscillons). We will establish the conditions 



for their existence, the reason for their longevity, and clarify their final collapse. Armed 
with a better understanding of their properties, we will also be able to suggest several 
situations where we believe oscillons can be of importance. 

The rest of this paper is organized as follows. In the next Section we will set up 
the general formalism and obtain the exact solution of the spherically-symmetric linear 
Klein-Gordon equation. As expected, in the linear case no oscillons appear, with bubbles 
quickly decaying away. We obtain the time-scale in which this decay occurs in order to 
later compare it to the case when nonlinearities are present. In Section 3 we present 
the numerical results that establish several of the key properties of oscillons for symmet- 
ric double-well potentials. Guided by these results, in Section 4 we present analytical 
arguments to explain why there is a minimum initial radius for bubbles to settle into 
the oscillon stage, why some oscillons live longer than others, and how oscillons finally 
disappear. In Section 5 we extend the numerical analysis of Section 3 to asymmetric 
double-well potentials, showing how the lifetime of oscillons is sensitive to the amount of 
asymmetry between the two minima. Here one must be careful to set the initial radius 
to be smaller than the critical radius, as bubbles with radii larger than critical will grow. 
As in the symmetric potential case there are no critical bubbles, we can say that we are 
studying the evolution of subcritical bubbles in symmetric and asymmetric potentials. 
Oscillons are thus a possible stage in the evolution of subcritical bubbles toward their 
demise. In Section 6 we discuss several possible situations in which these configurations 
will play an important role. Although we focus mainly on cosmological phase transi- 
tions, some of our arguments apply equally well to phase transitions in the laboratory. 
We conclude in Section 7 with a summary of our results and an outlook to future work. 
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II. PRELIMINARIES 



In this Section we introduce the notation and some definitions which will be useful 
later on. We also present the exact solution for the evolution of a "Gaussian-shaped" 
bubble (i.e. with 0(r, t = 0) ~ exp[— r 2 /R 2 ]) in the linear regime. 



A. General Formalism 

The action for a real scalar field in (3+l)-dimensions is 

"I 



S[*] 



d x 



;(^0)(^0)-Vs(A)(0) 



(1) 



where the subscripts S and A stand for symmetric (SDWP) and asymmetric (ADWP) 
double-well potentials, given respectively by, 

,2\ 2 



and, 



I7 , ,s TTl o ®0 m ,S ^ id 

\/ A (0) = — 2 - ^-0 3 + -0 4 . 



(2) 



(3) 



Note that the coupling constants A and a are dimensionless. A solution 0(x, t) to the 
equation of motion, 

dV(<f>) 



d 2 <p/dt z - V 2 



d<j) 



(4) 



has energy 



m = J 



d 3 x 



l(90/9t) 2 + i(V0) 2 + \/(0) 



(5) 



We will restrict our investigation to spherically-symmetric configurations. In this case 
it proves convenient to introduce dimensionless variables, p = rm, r = tm, and $ = ^4>- 
The nonlinear Klein-Gordon equation is, 
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^£ _ ^£ _ ^ _ I $ " $3 (SDWP) 

0r» dp 2 p dp ~ I _ $ + a$2 _ $3 (ADWP) 

where a = A _1/,2 aio. Note that for the SDWP the two minima are located at $o — — 1 
and $ + = 1. For the ADWP (with ct > 2), the minima are at $ — and $ + = 



1/2 



Requiring <3> + to be the global minimum implies a 2 > 9/2. For 
a 2 = 9/2 the two minima are degenerate with $ + = \fl. This value will be important 
later. 

We are interested in following the evolution of unstable spherically-symmetric config- 
urations of initial radius R , (from now on we call these initial configurations subcritical 
bubbles both for the SDWP and the ADWP) which can be thought of as being localized 
fluctuations about the global vacuum $ - Thus, we must measure the rate at which the 
initially localized energy is radiated away as the subcritical bubble relaxes to the global 
vacuum. This can be done by surrounding the initial configuration with a sphere of 
sufficiently large radius, R s ^> R , and measuring the flow of energy through the surface 
of the sphere. The evolution of subcritical bubbles is obtained by solving the nonlinear 
Klein-Gordon equation numerically. We define the subcritical bubble's kinetic, surface, 
and volume energies, respectively, by 

E k = 2ir [ Rs pV dp, E s = 2ir [** p 2 ($') 2 dp, E v = 4tt /"*" pV($) dp , (7) 

JO JO JO 

where a prime denotes derivative with respect to p. The bubble's total energy is thus, 

E b (r) = E k (r) + E s (t) + E v (t) . (8) 

In order to solve the nonlinear Klein-Gordon equation we will impose the following 
boundary conditions, 

$(p - oo, r) = $ , $'(0, r) = 0, $(p, 0) = . (9) 



The first condition guarantees that the bubble approaches the vacuum at <3> & t spatial 
infinity. The second condition imposes regularity at the origin, while the last condition 
states that the bubbles start their evolution at rest. These conditions must be supple- 
mented by the initial profile of the bubble. We will investigate both 'Gaussian' and 'tanh' 
bubbles which we write as 



$ c is the value of the field at the bubble's core, which we may or not take as being the 
other minimum of the potential, <J> + . If we do, the bubble can be interpreted as being a 
field configuration of initial linear size ~ 2R which interpolates between the two vacua. 
As we will see later, it is not necessary to set $ c = $ + in order to have subcritical bubbles 
relaxing into oscillons during their evolution. We can now move on to study the evolution 
of subcritical bubbles in the linear regime. 



As a first application of the above formalism, we will investigate the evolution of 
Gaussian bubbles in the linear regime. We choose as the linear potential, 



as it has a minimum at $ = — 1 with the same curvature as the SDWP. The Klein- 
Gordon equation has a trivial solution $(p, r) = —1. Separation of variables with a 
constant —k 2 allows us to write r) = — 1 + R(p)exp[±i\/ k 2 + 2r], with the radial 
function R(p) obeying, 




B. Bubble Evolution in the Linear Regime 



VL($) = ($ + 1) 2 



(12) 



R" + -R' + k 2 R = . 
P 



(13) 
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This equation has solutions which are linear combinations of sm kp , cos kp . Since cos kp is 



singular at the origin we write the general solution as 

sin kp 



$(p,r) 



poo I 

/ dkb(k)- 

Jo 



P 



cos(\/P + 2r) + a(k) sm(Vk 2 + 2r) 



The boundary conditions are, writing $ c = 2q — 1, q an arbitrary constant, 

$(p,r = 0) = 2g e- p2/ ^-l 
$(p — > oo, r) = -1 
$'(p = 0,r) = 
$(p,r = 0) = . 

Eq. [16| is trivially satisfied. Eq. [15] determines b(k), 

2goe-^= r d kb{k)^. 

Jo p 



Taking the sine transform we can write, 



b{k) = — ^Im 

7T 



I" dppe- p2/R °e ikp 
Jo 



The integral can be easily done and we obtain 

b(k) 



7T 



(14) 



(15) 
(16) 
(17) 
(18) 



(19) 



(20) 



sinkp 
P 



(21) 

vanishes as p — >• 0. Choosing 



Regularity at the origin is also guaranteed, as 
the bubble to start at rest implies that a(k) vanishes. Thus, the final solution satisfying 
all boundary conditions is, 



$ (P)T ) = _i + ^f r dkk 

V7T Jo 



sin kp 



cos(Vk 2 + 2r) 



(22) 



In Figure 1 we show a plot of this solution, for initial amplitude go — 1 an d radius 
R = 3. The bubble performs damped oscillations as it decays into $ — — 1- 



It is instructive to investigate the behavior of the bubble's core with time, 

$(p = 0, r) = -1 + ^3 r dk k 2 e - R ° k2/4 cos(v / FT2 r) . (23) 
V 71 " •'o 

The integral is dominated by small values of k, k ^ 2i?g . Thus, we can approximate 



(24) 



the argument of cos(v / ^ 2 + 2r) for i? ~ 2 and write, 

$(p = 0, r) = -1 + ^>Re [e^ /°° dfc /c 2 e" R ^ 2 /4 e W*ren 

y/TT L JO 

Performing the integral we obtain 

«C = 0, r ) = -1 + (i - £° mm cos (V2r + | tan- (ffi ) . (25) 

Thus, the amplitude at the core decays as r _3//2 , while the frequency becomes constant 
for r > _Rq/v^2- The envelope of the core's amplitude decays to 1/e of its initial value 
above <3> = — 1 in a time (units restored) 

ti/e ^ 1.18^m _1 . (26) 

In Fig. 2 we compare the above analytical approximation with the numerical solution of 
the Klein-Gordon equation (more details later). The excellent agreement gives support 
to the accuracy of the numerical methods used. 



III. EVOLUTION OF SUBCRITICAL BUBBLES IN SDWP: NUMERICAL 

RESULTS 

In this and the next Section we will restrict our analysis to bubbles in SDWP. Section 
5 will deal with bubble evolution for ADWPs. The equation of motion is 

$ _ $" _ = $ _ $3 ) (27) 
P 

with boundary conditions given by Eqs. [15| - [18]. (For tanh bubbles or any other initial 
bubble profile, just replace Eq. [15] by the appropriate choice.) This equation was solved 
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numerically using a finite difference scheme fourth order accurate in space and second 
order accurate in time. The radial dimension of the two-dimensional grid moved outwards 
with the speed of light in order to avoid any radiation from being reflected on the lattice 
boundary and thus interfering with the bubble's evolution within the grid (dynamically 
increasing simulation lattice). The alternative, a sufficiently long but static grid, is 
extremely time-consuming for long-lived oscillons. The resolution was typically set to 
Ap = 0.1 and Ar = 0.05; the total energy E b + 4vr J^p 2 [±$ 2 + |($') 2 + V{$)] dp is 
then conserved throughout the evolution to better than one part in 10 3 . Additionally, 
high resolution experiments (Ap = 0.01, Ar = 0.005) produced the same results, with an 
energy conservation of one part in 10 5 . We used R s = 10 for the SDWP, and R s = 15 for 
the ADWP. For reasons that will be made clear soon, we were not interested in bubbles 
of large initial radius. 

Figs. 3a and 3b show the energy of Gaussian and tanh bubbles for several initial radii 
-Ro- (More examples can be found in Ref. ||13|| .) In all the examples we took $ c = +1; the 
bubble interpolates between the two vacua. It is clear that the evolution of the bubbles 
is very sensitive to the value of Rq. An extensive investigation showed that Gaussian 
bubbles with R ~ 2.4 and R > 4.5 quickly disappear, radiating their initial energies to 
infinity. However, bubbles with 2.4 < R < 4.5, settle into a period of long-lived stability 
where practically no energy is radiated away. This stage in their evolution, which we call 
the oscillon stage, can have a duration approaching 10 3 — 10 4 m _1 , which is remarkably 
large compared to both short-lived bubbles and to the typical time-scales found for the 
linear potential. Although the range of values for R which fall into an oscillon stage is 
sensitive to the initial profile of the configuration, the same results are obtained for other 
initial bubbles, such as tanh bubbles. This supports our previous claim that oscillons can 
be viewed as a possible stage during the evolution of subcritical bubbles; after shedding 
a sufficient fraction of their initial energy, the subcritical bubbles enter the oscillon stage 
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which is characterized by an energy with a nearly constant value of ~ 43m/ A, regardless 
of their initial radius. 

The core value of the field, $(0,r), performs anharmonic oscillations as shown in 
Fig. 4a. In Fig. 4b we show a sequence of snapshots of an oscillon. Except when $(0, t) ~ 
—1, an oscillon configuration is very well approximated by a half-Gaussian. We define 
the effective radius of a localized field configuration by, 



C p 3 




dp 


C p 2 




dp 



Rett <S Rs ■ (2£ 



In Fig. 5 we show the evolution of this radius for several bubbles during the oscillon 
stage. The divergence at the end is spurious, signaling that no energy is left within the 
sphere (the denominator of Eq. ^). It is clear that the effective radius of an oscillon 
is approximately constant, with variations which are smaller than 20% about a mean of 
R osc ~ 2.8 — 3.0. This justifies the name given to these configurations: An oscillon is 
a localized, time-dependent field configuration with nearly constant radius and energy 
which is characterized by anharmonic oscillations of the field amplitude about the global 
vacuum. 

In order to stress the remarkable longevity of oscillons we show in Figs. 6a and 6b 
the lifetime as a function of initial radius and energy, respectively. For Gaussian bubbles 
the longest living oscillon, with n ~ 7.4 x 10 3 , comes from an initial bubble of radius 
Rq = 2.86. For tanh bubbles, the longest living oscillon comes from an initial bubble of 
radius Rq = 3.08, with lifetime T\ ~ 4.4 x 10 3 . In Fig. 7 we show the detailed dependence 
of lifetime as a function of radius for Gaussian bubbles about the peak at Rq = 2.86. 
(Lifetimes are accurate to within 5%.) 

So far we have restricted our investigation to bubbles that interpolate between the 
two vacua. This is it not a necessary condition for the existence of oscillons although, of 
course, it is sufficient. As long as the initial value of the field at the bubble's core probes 
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the nonlinearities of the potential and the initial radius is within the correct range (which 
varies with initial amplitude), oscillons can exist. We will give an analytical argument 
for this in the next Section. For now, we will just provide numerical evidence for this 
fact. In Fig. 8 we show a plot of lifetime for different core values. Clearly, no oscillon can 
develop if the initial energy is below the plateau energy. Also, we find that no oscillon 
develops if <3> c < $ in f, where <3> in f = — 1/V3 is the inflection point closest to $ - Thus, we 
arrive at the sufficient conditions for the existence of oscillons: i) the value of the field at 
the bubble's core must be above the inflection point, and ii) the initial bubble's energy 
must be above the plateau energy. Conditions i) and ii) fix the value of Rq for a given 
initial bubble to evolve into an oscillon. 

This concludes the presentation of our numerical results. In the next Section we will 
provide semi-analytical arguments to elucidate some of the properties of these configu- 
rations. 

IV. PROPERTIES OF OSCILLONS 

From the results of the previous Sections, it is clear that there are four main questions 
concerning the oscillons. First, why only bubbles with an initial radius above a certain 
value develop into oscillons, and how this value depends on the initial amplitude of the 
field at the bubble's core. Second, why certain oscillons live longer than others. Third, 
what is the mechanism responsible for the oscillon's final collapse. And finally, why above 
a maximum initial bubble radius no oscillons are possible. In this Section we address the 
first three questions. Work on the fourth question is in progress. 

In order to treat these questions analytically, we make use of the fact that inde- 
pendently of the initial bubble profile, an oscillon is very well approximated by a half- 
Gaussian. Even though the Klein-Gordon equation implies that $(p — > oo, r) ~ exp[— p], 
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the difference turns out to be sufficiently small in practice to justify our approximation. 
In a sense, the tail matters little to the dynamical properties of the configurations. The 
agreement of our analytical arguments with the numerical results should convince the 
reader of this fact. We model the oscillon by writing, for the SDWP, 



With this ansatz we have effectively reduced the field theory problem to two degrees of 
freedom, the amplitude at the core $ C ( T ) = 2q , (r) — 1, and the radius R{r). This problem 
is still quite complicated to treat analytically due to the nonlinear coupling between the 
two degrees of freedom. Further simplification is guided by the numerical investigation, 
which showed that the effective radius of the oscillon remains practically constant, with 
oscillations about its mean value of order 20% or less. Thus, as a first step, we will 
keep the radius constant, and treat only the amplitude at the core as an effective degree 
of freedom. Strong as it may seem, this simplification will allow us to extract several 
important results concerning the observed numerical behavior of these configurations, 
as we show in the next subsections. We are currently investigating the consequenses of 
keeping both degrees of freedom g(r) and R{r). 

The above model for the oscillon still misses one important ingredient; it does not 
include radiation of the bubble's energy to infinity. The justification for neglecting this 
lies in the fact that oscillons hardly radiate. By excluding radiation it is possible to 
analytically integrate the energy over the whole space. Using the definitions in Eq. |7| we 
obtain, for the kinetic, surface and volume energies, respectively, 



$(p,r) = 2g(r)exp -p 2 /R 2 {r) -1 



(29) 



7tv 2n 



R 3 q\ E t 



37TV27T 



Rq 2 , E v = ttV^R 3 q 2 




(30) 
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A. Existence of Oscillons: Lower Bound on the Initial Radius 



From Fig. 6a it is clear that there is a lower bound on the initial value of the bubble 
radius so that it relaxes into an oscillon during its collapse. Since from our previous 
discussion we know that oscillons are a product of the nonlinearities in the system, 
this result suggests that for small enough initial radii the nonlinearities are ineffective 
to trigger the resonant behavior responsible for the oscillon's longevity. That this is 
the case can be shown by studying the effective potential controlling the behavior of 
the amplitude q(r). Using the above ansatz with constant radius, the energy of the 
configuration E = + E s + E v can be written as the energy of a particle of unit mass 
with a potential V(q), 

E k + V{q), V(q) = (!+*)<? + , (31) 



16vrA 2" ^' V A) A A 

where A = ^J# 3 , B = ^pi?, C = ^R 3 , and D = >§R 3 follow from Eq. |3U|. The 
potential V has only one minimum at q — 0, ($ = $o = — 1, the global vacuum), about 
which the amplitude performs anharmonic oscillations. 

It is the energy localized within a small region surrounding the bubble that may (or 
not) feed the nonlinear growth of the modes ultimately responsible for the appearance of 
the oscillon during the collapse of the bubble. This lends further support to the above 
ansatz neglecting radiation. Thus, the equation of motion for the amplitude q(r) is, 

g = -2(l + !)g + 3^g 2 -4^g 3 . (32) 

Writing q{r) = q{r) +Sq(r), the linearized equation satisfied by the fluctuations Sq(r) is, 

5q = -u 2 (q,R)5q, oo 2 (q ,R) = 3^2f - 8^-q + (2 + A) , (33) 

where we have substituted the numerical values of the constants A, B, C, and D in the 
expression for the frequencies w 2 (q, R). Note that w 2 (q, R) is simply the curvature of the 

15 



potential dictating the dynamics of the amplitude q(r). As the bubble radiates its energy 
away, the configuration decays into the vacuum. However, for u 2 (q,R) < 0, fluctuations 
about q are unstable, driving the amplitude away from its vacuum value. These are 
the fluctuations which are mainly responsible for the appearance of the oscillon. In 
Fig. 9 we show a plot of the surface u 2 (q, R). It has one minimum at g m in ~ 0.77 (with 
location independent of R\), where its value is u 2 (q m i n ,R) ~ —0.514 + 3R~ 2 . Thus, 
only for R > R min ~ 2.42, u 2 < and fluctuations can grow. In other words, only for 
R > -R m i n — 2.42 are oscillons possible. This lower bound on the value of the radius 
agrees very well with our numerical results (see Fig. 6a). It is independent of the initial 
amplitude of the bubble. All bubbles with initial radius smaller than -R m i n will quickly 
collapse. (For core amplitudes above $ c = 1, it is possible to decrease the initial radius 
by about 20% or so and still obtain oscillons.) 



B. Collapse of Oscillons 

The above analysis can also provide information about the final decay of oscillons. 
For R > -Rmin, ui 2 (q, R) will be negative for amplitudes, 

4^3 



q- < q < q+ , q± 
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2R 2 



(34) 



Thus, for R > -R m m there is a minimum value for the amplitude at the core, shown in 
Fig. 10, $~(.R) = 2q_(R) — 1, below which the oscillon slips into the linear regime and 
quickly decays. This result can be understood as follows: As the bubble settles into the 
oscillon configuration with energy given by the plateau energy E ~ 43m/A and radius 
-R c fj ~ 2.8m -1 , there is a maximum value for the amplitude of the field at the core. 
This value is obtained from the formula for the static energy with its value fixed at the 
plateau value and with radius R ~ -R eff , and it is $ c ~ 0.2. The values of the field at 
the oscillon's core obtained numerically are always marginally within the allowed region 
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which gives uo 2 < 0; the oscillon survives while fluctuations are unstable. However, during 
the oscillon stage, energy is slowly being radiated away, and thus the amplitude at the 
core is slowly decreasing while the average value of the radius is slowly increasing (Fig. 5). 
From Fig. 10 and the argument above, below a certain value for the amplitude at the core 
the perturbations enter the linear regime and the oscillon decays. A comparison between 
$~(i? e fj) and the numerical values of the core's amplitude at the last oscillation is given 
in Table 1. In interpreting these results, we must keep in mind the crudeness of the 
analytical approximation used to obtain $~(i?). Even so, at least for the longest living 
oscillons, it is clear that during the last oscillation the amplitude falls below $~(i?). A 
more detailed analysis shows that the amplitude falls below Q~(R) during the last few 
oscillations, as the configuration starts to approach the linear regime responsible for its 
final demise. For completeness, in Fig. 11 we show a phase-space portrait of the evolution 
of ( &c(t) during the oscillon stage and its final collapse, for a bubble with initial radius 
Rq = 3.0. Clearly, the final spiraling into $o = ~ 1 ; typical of the linear regime, occurs 
as the maximum core amplitude (for $ = 0) falls roughly below $~. 



A question which is of great interest is the determination of the oscillon's lifetime as a 
function of the bubble's initial radius and core value. Although we were unable to obtain 
an analytical expression for the lifetime, we do understand why some oscillons live longer 
than others. Our argument is based on the virial theorem for spherically-symmetric 
scalar field configurations, which we derive next. 



Multiplication of the equation of motion, d 2 (p/dt 2 — d 2 (p/dr 2 — (2/r) d(p/dr = — 9V q^ 
(cf. Eq. |J), by 4irr 2 <j) and integrating over r gives, after integration by parts, 



C. Lifetime of Oscillons 




(35) 
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where we assumed that lim^oor 2 ^' = 0. The second term is easily recognized as twice 
the total surface energy. Performing a time averaging over one period, denoted by 

we get 

477- too r I f-oo fj\/ \ 

— r 2 dr dt(j>(j> + 2{E S ) + 4tt ( / r 2 0— dr ) = (37) 
T Jo Jt \Jq o<p J 

and after integrating by parts the time integral in the first term, (the boundary term (fxf> 
vanishes due to the integration over a period), 

-^J T dU7r J r 2 (f) 2 dr + 2(E s ) + 4tt M r 2 0— dr) = . (38) 

Identifying the first term as twice the time-averaged kinetic energy, we arrive at the virial 
theorem 

(E k ) = (E s )+2n^J\ 2 ( p^dr^ . (39) 

For the SDWP, in dimensionless variables, 

(E k ) = (E s ) + 2tt p 2 $ 2 ($ 2 - 1) dp^ . (40) 

As usual, the virial theorem holds as an equality only for strictly periodic systems. 
Numerical simulations of oscillons show, however, that the basic oscillation is overlaid by 
a long-wavelength modulation and other deviations from strict periodicity. It is hence of 
interest to analyze the "departure from virialization" , 

V(r) = (E k ) - (E s ) - 2n^j Q Rs p 2 $ 2 ($ 2 - 1) dp^j , (41) 

where now 

E k = 2ir [ S p 2 $ 2 dp, E s = 2tx p 2 ($') 2 dp (42) 

JO JO 
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and R s ~ 10 as before is an integration cut-off large enough to encompass the entire 
configuration. For a perfectly virialized configuration, V(t) = 0. In Fig. 12 we show the 
evolution of V for several Gaussian bubbles. When contrasted with Fig. 6a, it becomes 
clear that the longer the lifetime of the oscillon, the better virialized it is. This result is 
made more transparent by plotting the lifetime as a function of the maximum value of 
V for several radii, as shown in Fig. 13. Note also the symmetry about the longest-living 
oscillon, with R = 2.86. 

Using the virial relation and the numerical results, we can obtain a semi-analytical 
estimate for the optimal radius for an oscillon, that is, the one which is longest-lived. 
Although we perform the calculation for the SDWP, our methods can be easily generalized 
for any potential. With the ansatz for the Gaussian profile given in Eq. |29|, the time- 
averaged oscillon energy, and the departure from virialization, V, are, respectively, 

% = « 2 > + V2tf> - 4- i( 9 ">) R' + ^V>« , (43) 

and 

^ = + 2< 9 > - ^Itf> + - (,*>) *» - 3 -f(f)R . (44) 

Multiplying the expression for the time-averaged energy by 2 we can eliminate the cubic 
and quartic terms in the expression for V. Using that for the longest-lived oscillon V ~ 0, 
we obtain a cubic equation for the optimal radius, R mSlX , 

(^V> - + 2(g)) i?Lx + ^V>^ max - 2n- 3/2 (E) - . (45) 

To proceed, we further assume that q(r) is periodic, which is a good approximation for 
the longest-lived oscillon. Writing q{r) = q cos(u>r), with q an amplitude determined 
numerically (the reader should be careful to distinguish between this u and the one used 
in the linear perturbation analysis), the time-averaging can be performed and we finally 
obtain, 

19 



^ (3^ 2 - l) q 2 Rl ax + ^q 2 R max - 2tt- 3 / 2 (E) « . (46) 

The roots of this equation are determined once we know the values of the parameters 
(E), go, and uj. These can be obtained numerically using the remarkable independence of 
oscillons on initial conditions. We use (E) ~ 43, and uo = 2-k/T ~ 1.37. The maximum 
core amplitude, $ c , is roughly bounded by —0.1 ^ $ c ^ 0.2, which gives for q the range 
0.45 <• go ~ 0.6. With these parameters, we find that the equation has only one real root, 
bounded by 2.90 < -R ma x ~ 3.54. This range of values is in excellent agreement with the 
observed numerical range for the oscillon radius (see Fig. 5) providing strong support to 
our arguments. It also gives the correct range of initial values for the radius of bubbles 
which will relax into the longest-lived oscillons. Thus, the oscillon can be interpreted as 
the attractor field configuration which minimizes the departure from virialization. 

V. EVOLUTION OF SUBCRITICAL BUBBLES IN ADWP: NUMERICAL 

RESULTS 



It was first noted in Ref. |T3| that oscillons will also be present for nondegenerate 
potentials. Most of the analytical arguments above will also apply in this case. In 
particular, the minimum radius for subcritical bubbles to evolve into oscillons can also be 
obtained by the perturbation analysis presented in Section 4. A. The sufficient conditions 
for the existence of oscillons will still be the same, namely, that the initial energy be above 
the plateau energy, and that the initial core amplitude be above the inflection point of the 
potential. Of course, the plateau energy will depend on the degree of asymmetry of the 
potential. The important difference is that for ADWPs, the 0(3)-symmetric equations of 
motion admit static solutions known as bounces |16|]. These are the well-known critical 



bubbles of strong first order phase transitions, which specify the thermal barrier for the 



decay of metastable states, E cvlt |I7|; bubbles with radii larger than critical will grow, 
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converting the metastable phase into the stable phase with lower free-energy density. 
Thus, when discussing oscillons in the context of ADWPs, we must make sure that the 
initial configurations have radii smaller than the critical bubble radius, R cr it, as well as 
energies smaller than the decay barrier. The initial bubble energy is bounded by the 
plateau energy from below and the decay barrier from above. 

In order to see the effects of the asymmetry on the properties of the oscillons, we start 
by showing the results for the degenerate case, obtained by setting a = 3/ a/2 in Eq. |[ 
Recall that in this case the minima are at $o — and $ + = a/2- In Fig. 14 we show 
the lifetime of oscillons as a function of initial radius for several core amplitudes. Note 
that the lifetimes are larger than for the SDWP (Fig. 6a). This is simply due to the fact 
that for a = 3/ a/2 the ADWP is shallower and narrower than the SDWP, softening the 
surface energy of the initial bubbles. In Fig. 15a we show the lifetimes vs. radii of initial 
Gaussian bubbles leading to long-lived oscillons for different values of a. For reference 
we also show the values of the critical radii. The perturbative analysis of Section 4. A can 
easily be adapted to this ADWP case, yielding an expression for the frequencies u> 2 (q, R) 
of small fluctuations (analogous to Eq. [33]) , 



The minimum of this surface (for fixed a) is once again independent of R with 
q min ~ 0.51a, hence oscillons are possible only for _R > -Rmin — (3/(0. 28a 2 — l)) 1 ^ 2 . 
With the values a = 3/a/2, 2.16 and 2.23 we then obtain R > 3.39, 3.12 and 2.76 
respectively, results which compare favorably with the numerical simulation values of 
Ro > 3.2, 3.1 and 2.9, respectively. In Fig. 15b we show lifetime vs. initial bubble energy 
for different values of a. For reference we give the values of the plateau energy and of the 
decay barrier. Note that as the asymmetry is increased, the lifetimes of the oscillons also 
increase, almost by a factor of two between the nearly degenerate a = 2.16 and the more 




(47) 
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asymmetric a = 2.23, while the ratio between the critical bubble radius, -R cr it, and the 
longest-lived oscillon, -R max , varies from -R cr it/-R m ax — 5 for a = 2.16 to -R cr it/-Rmax — 2 
for a = 2.23. As the asymmetry is increased, the oscillons approach more and more the 
critical bubble, explaining their increased longevity. 

VI. OSCILLONS IN ACTION: POSSIBLE APPLICATIONS 

In this Section we will present a few situations in which we expect oscillons to be 
relevant. As we will argue below, their remarkable longevity makes them specially in- 
teresting in the context of phase transitions; if thermally nucleated, their presence can 
affect the dynamics of the transition in several ways. It is not our intention here to give 
a detailed treatment of the role of oscillons on the dynamics of phase transitions, but 
simply to stress the interesting physics that can emerge due to these configurations. 

As we have seen, a typical range of lifetimes is between ti = 10 3 — 10 4 m" 1 in units of 
the mass m introduced in Eqs. ^| and |3|. This is much longer than that of the solution to 
the spherically symmetric linear Klein-Gordon equation ~ 5m" 1 . The expansion rate of 
the Universe in a radiation-dominated regime can be written in terms of the background 
temperature T as H 2 oc T 4 /m P1 , where H is Hubble's parameter. Thus, the expansion 
time-scale is tu ~ H~ x oc (mpi/T)T _1 . Typically, the symmetry breaking temperature 
T c can be written in terms of the mass scale m of the theory as T c ~ m/ \f\. Thus, the 
expansion time-scale at T c is tu ~ A(mpi/m)m _1 . The ratio between the oscillon lifetime 
(taking ti = 10 4 m _1 ) and the expansion time-scale is then, ti/tu ~ A _1 10 4 (m/mpi). 
From this we see that for masses of order the GUT scale the lifetime of the oscillons is 
comparable (or larger, for weak coupling!) to the age of the Universe at that scale, an 
intriguing possibility. In such a scenario these unstable field configurations could have 
a dramatic effect on the dynamics of any phase transition. For example, during a first 
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order phase transition those subcritical bubbles which go on to form oscillons could last 
long enough to become critical bubbles as the Universe cooled. In this way we would 
have a method of completing the transition quicker. Another possibility is that oscillons 
act as seeds or nucleation sites for the critical bubbles. The combination of these two 
effects will increase the production rate of critical bubbles, a feature which may well have 
useful consequences for the old inflationary universe scenario. That model failed partly 
because the production rate of critical sized bubbles and larger could not keep pace with 
the exponential expansion of the Universe. If the bubble nucleation rate were increased, 
then this problem may well be overcome. 

For oscillons to be relevant cosmologically, not only must they survive for long enough, 
they must be thermally produced in large enough numbers since they are unstable and 
eventually decay. A naive estimate of this rate is that the number density of oscillons of 
size R produced at temperature T due to thermal fluctuations is 

n(R,T)^T 3 e- F ^ T , 

where F(R) is the free-energy of the configuration of radius R and is given by F(R) = 
E s + E v in Eq. [30[ Comparing F(R osc ) with F(R cr ) gives an indication of the fraction 
of bubbles which are oscillons as opposed to critical at any given temperature T. In 
fact we can see quite easily that although the oscillons are unstable they are produced 
in much greater abundance, as their free-energy barrier is typically smaller than that 
for critical bubbles. To be sure of this we require that their thermal nucleation rate be 
considerably larger than the expansion rate of the Universe, i.e. T t ^(R,T)/ H ^> 1. Since 
H oc T 2 /mpi, it becomes a straightforward comparison. For a GUT scale transition, 
say with T c ~ 10 15 GeV, this condition implies F osc /T c < 10, which is not difficult to 
satisfy for sufficiently weak transitions, such as the SU(5) Coleman- Weinberg model, as 
shown in Ref . . More speculatively, even for lower energies, oscillons may be potentially 
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relevant. Although in these cases their lifetime becomes small compared to the expansion 
time-scale, if they are produced in large enough numbers these sub-critical bubbles could 
still affect the dynamics of the transition. The presence of oscillons will substantially 
increase the equilibrium number- density of sub-critical bubbles of the broken phase; as 
they last longer, their depletion rate by shrinking is decreased, and hence the net volume 
in the broken phase increases. (In Ref. |ll|], it was assumed that all subcritical bubbles 



were roughly of a correlation volume and disappeared in a time ~ £,{T), where 
is the correlation length.) They effectively make the transition weaker than what one 
would predict from the effective potential. Also, the collision of oscillons with expanding 
critical bubbles will possibly cause instabilities on the bubble wall, implying that the 
assumption of spherical evolution of the walls may be incorrect. For the case of the 
electroweak transition, Gleiser and Kolb |9| have shown that the condition on subcritical 
bubbles can be written as F(R)/T < 34. This in turn imposes a constraint on the mass 
of the associated Higgs, which for the one-loop potential turns out to be mHi ggs > 88GeV 
(see JIB] for details). It seems to be the case that oscillons will have an important effect 
on the dynamics of sufficiently weak first order phase transitions. 

In this paper we have been investigating the existence of oscillons for both first and 
second order phase transitions. A number of issues arise common to both cases which 
require further study. The first concerns the coupling of oscillons to other forms of 
matter, whether they be other scalar fields, gauge fields or fermions. We have regarded 
the oscillons as emerging from an effective theory in which the fields to which it is coupled 
have been integrated out (a procedure we would hope is valid for low enough energies), 
leaving an effective potential for the scalar field. Ideally we would like to consider the 
full theory and solve for all the fields without integrating out the massive ones. It could 
be that one of these fields leads to an instability in the $ field which causes the oscillon 
to decay faster than we have estimated. On the other hand, coupling the oscillon to 
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a charged field may enhance its lifetime, as in the case of nontopological solitons [|J. 
A second issue concerns the coupling of oscillons to hot plasmas, as would be the case 
during thermal phase transitions. The plasma would act both as a viscous medium and 
as an enhancer of fluctuations, presumably affecting the lifetime of the oscillons. We are 
currently investigating both issues. 

The discussion in this Section has concentrated on early Universe aspects of oscillons. 
Since they are field theories we should expect them to be seen at laboratory energies as 
well. This may not be so easy to do in practice, but there are many examples of phase 
transitions in liquid crystals and Ising-like systems which produce nonlinear field theory 
objects such as topological defects ]n|. Also, solutions to the nonlinear Schrodinger 
equation have been known to be of importance in several contexts, including the propa- 
gation of information in optical fibres |2(| . It is reasonable to expect that oscillons will be 
present in the non-relativistic limit, thus being possible solutions to the time-dependent 
nonlinear Schrodinger equation as well. What would be required for oscillons is a distinct 
signature. It could be that as the energy of the system reaches its plateau during the 
oscillon stage, the material has a particular refractive index and thus could be detected 
in scattering experiments. 

VII. CONCLUSIONS AND OUTLOOK 

In this paper we have presented the results of a detailed investigation of the properties 
of oscillon configurations and explained, where possible, the physics behind their interest- 
ing dynamics. The fact that they exist in both first and second order phase transitions 
makes them of particular interest. They are localized, non-singular, time-dependent, 
spherically-symmetric solutions of nonlinear scalar field theories, which are unstable but 
extremely long-lived, with lifetimes of order 10 3 — 10 4 m -1 , where m is the mass of the 
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scalar field. They naturally appear during the collapse of spherically symmetric field con- 
figurations. We have obtained the conditions required for their existence, namely that 
the initial energy needs to be above a plateau energy and the initial amplitude of the field 
needs to be above the inflection point on the potential in order to probe the nonlinearities 
of the theory (but does not need to be at the true minimum of the potential). 

Of the many intriguing aspects of these configurations, some that stand out include 
the fact that they exist only for a given range of initial radii and core amplitudes. The 
lower value of the radii can be explained by perturbation theory. It corresponds to the 
minimum radius beyond which the field probes the nonlinearity of the potential. Explain- 
ing the upper bound for the initial radius of the field profile is not so straightforward and 
we are currently investigating this. It could well be that since larger bubbles have larger 
initial energies, during their collapse higher nonspherical modes are excited, triggering 
the rapid growth of instabilities responsible for the bubble's collapse before it can settle 
into the oscillon stage. Another remarkable feature is that the plateau energy of the 
oscillon is practically independent of the initial radius. We have interpreted this fact by 
showing that the oscillon can be thought of as the attractor field configuration which 
minimizes the departure from virialization. 

There is much that remains to be investigated. One concern is that we only investi- 
gated stability to radial perturbations. We really need to investigate how nonspherical 
perturbations affect the spherically symmetric solutions. One possibility is that they 
will tend to make the oscillons collapse into a pancake configuration, and hence decay 
more quickly than in the spherical case, although we believe this will only be the case 
for bubbles with large initial radii. We may also think of higher nonspherical modes as 
excited states of the "ground-state" I = resonance studied here. It is thus possible 
that oscillons may appear in higher energy configurations, which may decay either to the 
ground-state oscillon or just into scalar radiation. Finally, a more detailed study of the 
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coupling of these objects to other matter fields and hot plasmas is required in order to 
investigate how they affect the dynamics of phase transitions and how their own decay 
is affected by these couplings. It is clear though that they are of interest cosmologically. 
We are currently analyzing the consequences of oscillons if they were to be formed at the 
electroweak scale |18| . 
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